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TECHNICAL NOTE 2892 


A RAPID METHOD FOR ESTIMATING THE SEPARATION POINT OF A 
COMPRESSIBLE LAMH*AR BOUNDARY LAYER 
By Laurence K. Loftin, Jr., and Homer B. Wilson, Jr. 


SUMMARY 


A method has "been developed for rapidly estimating the separation 
point of a laminar boundary layer in a compressible flow. The method 
consists of an extension of Von Doenhof f ' s simplified solution for the 
incompressible case (NACA TN 671) and makes use of a set of transforms 
derived by Stewartson (Proc. Roy. Soc., 19^9) which permit, under cer- 
tain assumptions, the expression of compressible laminar layers in termB 
of equivalent incompressible laminar layers. The method developed is 
generally applicable to any two-dimensional flow in which the classical 
boundary- layer assumptions are satisfied. The dependence of the method 
upon the boundary- layer assumptions, of course, means that it should 
not be applied to determine whether the pressure rise through a shock 
wave causes separation. 

Calculations of the laminar separation point for a wide range of 
Mach number and velocity gradient indicate that, for all velocity gra- 
dients, the amount of velocity recovery possible before laminar separa- 
tion occurs decreases as the Mach number increases. 


INTRODUCTION 


The determination of the position at which the laminar boundary 
layer separates has been the subject of much theoretical investigation. 

For the case of incompressible flow, such methods as those of Howarth 
(ref. 1), Von Karman and Millikan (ref. 2), and Von Doenhoff's simpli- 
fied adaptation of the Von Karman-Mi lli kan theory (ref. 3) may be used 
for determining the laminar separation point on bodies with velocity 
distributions of arbitrary shape. 

For the case of compressible flows, several investigators have been 
interested in determining the basic nature of the effect of Mach number 
on the laminar separation point. Such methods as those of Howarth (ref. 4) 
and Cope and Hartree (ref. 5) provide means for determining the laminar 
separation point in a compressible flow. These methods are rather 
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cumbersome, however, and are not readily adaptable for rapidly estimating 
the effect of Mach number and velocity gradient on 1 am-t nar separation. 
Stewartson (ref. 6) has recently shown that a vast simplification in 
the theory of compressible laminar boundary layers is possible if cer- 
tain rather reasonable simplifying assumptions are made. On the basis 
of these assumptions, he was able to develop transforms which express 
the compressible 1 a mi n a r boundary layer in terms of an equivalent 
incompressible 1 ami nar boundary layer. 

In the present paper the method of Von Doenhoff has been combined 
with Stewartson' s transforms to provide a means for rapidly estimating 
the separation point of a laminar boundary layer in a compressible flow. 
The method developed is generally applicable to any two-dimensional flow 
in which the classical boundary-layer assumptions are satisfied. With 
the -use of the method developed, the effect of Mach number on the laminar 
separation point has been calculated for a wide range of Mach number and 
velocity gradient. The results of these calculations are presented. 


SYMBOLS 
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reference length 

abscissa 

ordinate 

equivalent flat-plate length 
speed of sound 

velocity just outside boundary layer 
velocity inside boundary layer 
Mach number 

L dU 

nondimensional velocity gradient, 

u o dx 

specific heat at constant pressure 
specific heat at constant volume 
ratio of specific heats, Cp/c v 


viscosity 
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p mass density 

V kinematic viscosity 

k coefficient of heat conduction 

Subscripts: 

0 reference condition taken at point of maximum velocity 

1 any point along body just outside boundary layer 

c compressible 

i incompressible 


ANALYSIS 

Stewartson's Transforms 


Stewartson (ref. 6) showed that the equations governing the behavior 
of two-dimensional laminar boundary layers in a compressible flow are 
identical to those governing their behavior in an incompressible flow if 
the variables in the two planes a re related by the following transforma- 
tions: 


y± = 


ry c 

<Wvo 0 



(1) 


37-1 



(3) 
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where x<>, y c , U c and xq, yq, Uq are the coordinates and the velocity 
just outside the boundary layer in the compressible and incompressible 
planes , respectively. The subscript o refers to some convenient refer- 
ence condition and the subscript 1 refers to conditions at a variable 
point just outside the boundary layer. The detailed derivation of the 
incompressible-boundary-layer equation from the compressible equations 
with the use of transforms (l) to (3) is treated in reference 6. 

The transformation of the compressible-boundary- layer equation to 
an equivalent incompressible-boundary- layer equation with the use of 
transforms (l) to (3) depends upon the following assumptions: The 

boundary is thermally insulated, the viscosity is proportional to the 
absolute temperature, and the Prandtl number Cpd/k of the fluid is 
unity. By assuming the boundary thermally insulated, the basic problem 
of motion of the compressible boundary layer is isolated from extraneous 
problems associated with heat transfer through the boundary surface. 

The assumption that the Prandtl number is unity means that the stagna- 
tion temperature is reached at the wall. The Prandtl number for air, 
however, is actually about 0.7T5 (ref. 4); therefore, the temperature 
near the wall is lower than that predicted. Since the viscosity 
increases with the temperature, the assumption that the Prandtl number 
is unity would indicate that the effects of viscosity are overestimated. 
With regard to the remaining assumption that the viscosity is proportional 
to the absolute temperature, experiments have shown that the viscosity 
for air varies as the eight-ninths power of the absolute temperature 
between 90° and 300° Kelvin (ref. A). Thus, this assumption also results 
in an overestimation of the effects of viscosity. 


Von Doenhoff 's Method 

Since the transforms (l) to (3) permit the expression of the com- 
pressible laminar boundary layer in terms of an equivalent incompressible 
laminar boundary layer, it is apparent that all methods for calculating 
the laminar separation point in incompressible flow are equally applicable 
in the case of compressible flow. The laminar separation point is 

8u 

defined as the point along the surface at which — =0 at y = 0 . Of 

oy 

all methods available for calculating the laminar separation point in 
incompressible flow, that due to Von Doenhoff (ref. 3) is perhaps the 
simplest over-all method to apply and, is used herein as a basis for 
developing a method for rapidly estimating the laminar separation point 
in a compressible flow. 


With the use of the theory developed by Von Kalman and Millikan for 
confuting the separation point of a laminar boundary layer. Von Doenhoff 
calculated the laminar separation point for a series of velocity dis- 
tributions which consisted of a region of uniform velocity (flat-plate 
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flow) followed by a region of linearly decreasing velocity. These 
calculations were generalized to the case of bodies with velocity dis- 
tributions of arbitrary shape on the basis of the fundamental assump- 
tion that the shape of the boundary- layer velocity distribution at the 
point along a body at which the outside velocity is a maximum , that is, 
at the point of application of the adverse velocity gradient, is very 
nearly the same as that of a Blasius flat-plate distribution. Thus, 
the condition of the boundary layer at the point of application of the 
adverse velocity gradient can be represented by the flow over an equiva- 
lent length of flat plate with a uniform velocity equal to the maximum 
velocity. The assumption was also made that, for purposes of estimating 
the laminar separation point, any velocity gradient likely to be encoun- 
tered could be approximated by a straight line. Von Doenhoff presented 
the results of his calculations in the form of the decrement in velocity 
necessary to cause separation AUi]_/Ui 0 as a function of the nondimen- 

L 1 dU ± 

sional velocity gradient Fj_ = — where Lq is the equivalent 

u io ax ± 

flat-plate length, Ui 0 is the maximum velocity, and dUi^/dxq is the 
velocity gradient. The relationship between AUi^/Ui Q and F^ given 

in reference 3 is presented herein as figure 1. The equivalent length 
of flat plate corresponding to the flow at the point of application of 
the adverse velocity gradient can be found with the use of the following 
equation (ref. 3, eq. (l) ) : 



(M 


where Uii/di 0 is the velocity ratio at any point on the body and the 

integration is carried from the beginning of the flow to the point of 
application of the adverse gradient. Equation (4) was obtained from 
an integration of the Von Karman momentum relation by making the 
assumption that the boundary- layer velocity distribution at every point 
along the body in the region of accelerating flow is the same as a 
Blasius flat-plate distribution. Calculations of the laminar separa- 
tion point on the NACA 0012 airfoil at zero lift made by the use of 
Von Doenhoff 1 s rapid method and calculations made by the use of the more 
elaborate Von Karman -Millikan theory are shown in reference 3 "to agree 
very closely. 
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Application of Stewartson's Transforms to Von Doenhoff 's Method 


In order to find the position on a body at which laminar separation 
occurs, the distribution of velocity outside the boundary layer and the 
Mach n umb er at some typical point must be known. With this information, 
the 1 nmi nnr separation point can be calculated through the direct appli- 
cation of transforms (2) and (3) and Von Doenhoff 's relation given in 
figure 1. Although relatively easy to apply, this method does involve 
a rather tedious point-by-point transformation from the compressible to 
the equivalent incompressible plane. 


For this reason, a more rapid procedure has been developed for 
estimating the separation point of a compressible laminar boundary 
layer. This procedure consists essentially in the direct application 
of Von Do enho ff ' « method to the compressible velocity distribution, the 
effect of Mach number being accounted for by a simple multiplying fac- 
tor applied to the measured velocity gradient in the compressible plane. 


Consider the velocity gradient dU c /dx c in the compressible plane. 

37-1 


From equation (2) it can be seen that dxi 


- ffiY" 1 

\ a o/ 


dx c and equation (3) 


states that Up, = — U c . Consequently, the compressible and incompres- 

£l-j^ -L 

sible velocity gradients are related as follows: 


or 
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If the reference velocity of sound ao is taken at the point of maximum 
velocity U Cq , that is, at the point of application of the adverse veloc- 
ity gradient, then = Ui 0 a n d 



By the use of the one-dimensional energy equation, the terms a^/a,-, and 
can be written in terms of the Mach number Mq and the 
velocity ratio U c ^/u Co : 




and 



An examination of equations (5) to ( 7 ) indicates that a linear velocity 
distribution in the compressible plane transforms into a nonlinear dis- 
tribution in the equivalent incompressible plane. The basic assumption 
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is made, however, that the determination of the appro xim ate laminar 
separation point in the equivalent incompressible plane can he made 
from a knowledge of the transformed adverse velocity gradient at the 
point of maximum velocity. This assumption presupposes that there is 
a discontinuity in the compressible velocity distribution at the point 
of application of the adverse gradient or, at least, that the velocity 
distribution can be approximated in this manner. On the basis of this 

a l 

assumption, — = 1.0 and 



Substitution of equation (8) into equation (5) and the value 1.0 for 
the terms U Cl /Uc 0 and &o/a.-^ yields the following relation: 



(9) 


Thus, the value of the velocity gradient at the point of maximum veloc- 
ity in the incompressible plane is obtained by multiplying the measured 
value of the velocity gradient in the compressible plane by the fac- 
tor 1 + 7 — Mo 2 , where it should be remembered that Mq is the Mach 
number at the point of application of the adverse gradient in the com- 
pressible plane. The factor 1 + 7 ~ 1 Mq 2 is plotted as a function 

u 

of Mach number in figure 2. 


The velocity gradient in the equivalent incompressible plane (eq. ( 9 )) 
must be made nondimens ional in terms of the equivalent flat-plate length 
in the incompressible plane if it is to be used with the relation given 
in figure 1 for finding the velocity recovery before laminar separation 
occurs in an incompressible flow. The equivalent length of flat plate in 
the incompressible plane can be determined from the relation 
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Li 




8.17 

dXjL 


( 10 ) 


or, in terms of variables in the compressible plane, from the relation 


8.17 8.17 


L = r C ° /5o\ * (%) 

\°<J W ^ 


37-1 

,7-1 


dx r 


(11) 


For air, y = 1.4; therefore, 


37-1 

7-1 


8.0 and equation (ll) reduces to 


Li = 



( 12 ) 


Equation (12) can be used directly with the compressible velocity distribu- 
tion to determine Lp. In many cases, however, it would appear that the 

/M 0 ’ 17 

value of I — would be so near unity as to have a negligible effect 

\ a l/ 

on Li- In this cane, the equivalent flat-plate length L^ can be 
determined from the equation 


Li = 




8.17 
) tac 


(13) 


The complete expression for the nondimens ional velocity gradient in the 
equivalent incompressible plane in terms of variables in the compressible 
plane is therefore 
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In the determination 

course, he expressed 
plane. 

The determination of the laminar separation point corresponding to 
a known velocity distribution at a given Mach number involves the fol- 
lowing steps: 

(a) Graphical measurement of the velocity gradient dU c /dxc in 

the compressible plane. This process may require fairing of the given 
velocity distribution in such a way that the adverse velocity gradient 
begins at a distinct point and is approximated by a straight line. 

(b) Calculation of the equivalent flat-plate length by the use of 
equation (13). 


dUi 


dxi 


i - A + 111 ^2L 

2 * M a dx c 


-(> 


(HO 


of Lp and 


^U Cr 


/dx c the value of x c can, of 


in terms of a reference length in the compressible 


(c) Determination of the factor 



from figure 2. 


The nondimens ional velocity gradient in the equivalent incompressible 
plane may then be determined by the use of equation (lk) . The velocity 
ratio corresponding to laminar separation in the incompressible plane 
may be found from the curve of figure 1 which expresses the velocity 
decrement necessary to cause laminar separation in terms of the non- 
dimensional velocity gradient. The velocity ratio corresponding to 
laminar separation in the compressible plane is found with the use of 
U ± , ^ U c x 

the relation — — = . The 

Ui 0 a l Uc 0 

Ucp/Uc 0 for different Mach numbers is plotted in figure 3 * 


relationship between Ui^/Ui 0 and 


In some cases, considerable fairing of the given velocity distribu- 
tion may be required in order to obtain what appears to be a reasonable 



dx c . Under such circumstances, the position of the 


measure of 


calculated laminar separation point with respect to the given velocity 
distribution may indicate that a somewhat different fairing of the 
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velocity distribution would give a more realistic approximation of that 
portion of the velocity gradient along which laminar separation occurs . 
A second calculation of the laminar separation point may then be 
indicated . 


Comparison With Results of Previous Analyses 

The effect of Mach number on the position at which laminar separa- 
tion occurs ip a boundary- layer flow progressing against a linearly 
decreasing velocity from the leading edge has been investigated by 
Howarth (ref. if) who employed a Von Rarman-Pohlhausen type of analysis, 
and by Stewartson (ref. 6) who applied his transforms to the incompres- 
sible solution developed by Howarth in 1938 (ref. l). This type of 
velocity distribution has a value of F c equal to zero since the 
equivalent length of flat plate is zero. In order to provide a com- 
parison between the method of the present investigation and those used 
in references if and 6 , the laminar separation point has been calculated 
by the method presented herein for the case of F c equal to zero for 
a range of Mach number. The results of these confutations are presented 
in figure b together with the results obtained by Howarth and Stewartson. 
Examination of figure b indicates that considerable difference exists 
at all Mach numbers between the predictions of the three methods. For 
zero Mach number, the value of the decrement in velocity ratio of 0.12 
given by Howarth in reference 1 is probably the most exact. The magni- 
tude of the differences between the three methods, however, appears to 
be roughly the same for all Mach numbers. This result can be seen more 
clearly in figure 5 in which the difference in AU Cl /u Co at Mach num- 
ber zero and at some arbitrary Mach number is plotted against Mach num- 
ber. The trends shown in figure 5 indicate that the predicted effect 
of Mach number on the laminar separation point as determined by the 
method of the present investigation is in substantial agreement with 
the results obtained by Stewartson and Howarth for the case of linearly 
decreasing velocity from the leading edge. 


EFFECT OF MACH NUMBER ON THE LAMINAR SEPARATION POINT 


In order to show the effect of Mach number on the position of 
laminar separation for a series of values of F c , the nondimensional 
compressible velocity gradient, calculations have been made for Mach 
numbers varying from 0 to 10 and, for nondimensional con^pressible veloc- 
ity gradients F c varying from 0 to -0.12. These calculations are for 
velocity distributions represented by a region of uniform velocity equal 
to the maximum velocity followed by a region of uniformly decreasing 
velocity. The results of these calculations are presented in figure 6 
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and show that the velocity recovery obtainable before laminar separation 
decreases as the Mach number increases ' for all nondimens ional velocity 
gradients. For a given velocity gradient, therefore, the lami nar separa- 
tion point moves forward as the Mach number increases. Figure 6 also 
shows that the shape of the curves of velocity recovery against Mach 
number is not influenced to any large extent by the value of the velocity 
gradient. 

It is important to point out that the individual curves presented 
in figure 6 show the effect on the .1 am in ar separation point of increasing 
the Mach number for fixed values of the velocity gradient. The problem 
of most practical interest, however, is that of a body having fixed 
geometry. In this case, the velocity gradients vary as the Mach number 
varies. Thus, when the laminar separation point on a given body for 
different Mach numbers is calculated, the value of the compressible 
velocity gradient F c varies with Mach number. In general, adverse 
velocity gradients become steeper as the Mach number increases so that, 
for a given body, the forward movement of the laminar separation point 
with increasing Mach number would be greater than is indicated in fig- 
ure 6 for a given velocity gradient. 

A word of caution should be added with regard to the type of flow 
field to which the method developed herein should be applied. It will 
be recalled that the method developed is based on the classical boundary- 
layer assumptions . These assumptions state that the rate of change of 
velocity in the x-direction must be small, with respect to itB rate of 
change in the y-direction. The method obviously cannot be applied in 
the vicinity of a discontinuity in velocity (pressure) such as that 
associated with the presence of a shock wave. 


CONCLUDING REMARKS 


A method has been developed for rapidly estimating the separation 
point of a laminar boundary layer in a compressible flow. The method 
consists of an extension of Von Doenhoff 's simplified solution for the 
incompressible case (MCA TN 671 ) and makes use of a set of transforms 
derived by Stewartson (Proc. Roy. Soc., 19^9) which permit, under cer- 
tain assumptions, the expression of compressible laminar layers in terms 
of equivalent incompressible laminar layers. The method developed is 
generally applicable to any two-dimensional flow in which the classical 
boundary- layer assumptions are satisfied. The dependence of the method 
upon the boundary- layer assumptions, 'of course, means that it should 
not be applied to determine whether the pressure rise through a shock 
wave causes separation. 
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Calculations of the laminar separation point for a wide range of 
Mach number and velocity gradient indicate that, for all velocity 
gradients, the amount of velocity recovery possible before laminar 
separation occurs decreases as the Mach number increases. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., November 19, 1952. 
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Figure 1. - Incompressible velocity decrement £SJ±^J\J± o obtainable before 

laminar separation occurs as a function of the velocity gradient Fi 
(ref. 3). 







Figure 4. - Comparison of the results of the present method vith those 
of Howarth and Stewartson for the case when F - 0. 
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Figure 6.- Compressible 

before separation as 
gradients F c . 


velocity decrement AU c ^jU c obtainable 
a function of Mach number for various velocity 
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